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Climate change will lead to loss of range for many species, and thus to loss of genetic diversity crucial for 
their long-term persistence. We analysed range-wide genetic diversity (amplified fragment length poly- 
morphisms) in 9581 samples from 1200 populations of 27 northern plant species, to assess genetic 
consequences of range reduction and potential association with species traits. We used species distri- 
bution modelling (SDM, eight techniques, two global circulation models and two emission scenarios) 
to predict loss of range and genetic diversity by 2080. Loss of genetic diversity varied considerably 
among species, and this variation could be explained by dispersal adaptation (up to 57%) and by genetic 
differentiation among populations (F ST ; up to 61%). Herbs lacking adaptations for long-distance disper- 
sal were estimated to lose genetic diversity at higher rate than dwarf shrubs adapted to long-distance 
dispersal. The expected range reduction in these 27 northern species was larger than reported for tem- 
perate plants, and all were predicted to lose genetic diversity according to at least one scenario. SDM 
combined with F S t estimates and/or with species trait information thus allows the prediction of species' 
vulnerability to climate change, aiding rational prioritization of conservation efforts. 
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1. INTRODUCTION 

When addressing impacts of climate change on biological 
diversity, most studies treat a species as a unit and thus 
ignore intraspecinc genetic variation [1,2]. Maintaining 
genetic diversity within a species is crucial for its ability 
to adapt both in the short-term and long-term survival 
[3-5]. Species may respond to climate change by local 
adaptation [6,7], range shift [8,9], range reduction [10] 
or a combination of these [6]. While range shift may 
alter the genetic diversity within species [11,12], range 
reduction is most likely to cause loss of genetic diversity 
[13,14] and may therefore severely limit the species' 
ability to adapt to a changing climate [4] . 

Most species adapted to cold environments are 
expected to suffer range reduction following climate 
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warming [7,10]. In animals, demographically challenged 
populations showed 22-26% reduction in genetic diver- 
sity relative to healthy populations [15], and two of 
three studied lizard species and one cold-adapted aquatic 
mayfly were expected to experience genetic depaupera- 
tion owing to climate change [16,17]. By contrast, in 
the only plant species studied to date, the wind-dispersed 
dwarf shrub Salix herbacea, a loss of 50 per cent of its 
European range was estimated to cause a loss of only 5 
per cent of its genetic diversity because of its high disper- 
sal ability and history of broad-fronted postglacial 
colonization [18]. Thus, species traits that influence the 
distribution of genetic diversity within and among popu- 
lations [19,20] may also determine the susceptibility of 
species to genetic diversity loss when their ranges are 
reduced, and their consequent vulnerability to loss of 
evolutionary potential [3] . 

Here, we use published and new amplified fragment 
length polymorphism (AFLP) data combined with a ran- 
domization procedure to estimate loss of genetic diversity 
under increasing loss of range for 27 northern plant 
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species, and explore the effect of different species traits. 
We also explore whether genetic differentiation among 
populations (F S t)> which is available for thousands of 
species, can be used as a proxy for predicting the genetic 
vulnerability of species to climate change. Finally, we use 
an ensemble of species distribution models (SDMs) to 
estimate the expected range reduction and associated 
loss of genetic diversity under two different global circula- 
tion models (GCMs) and two emission scenarios. Our 
results show that the expected genetic consequences of 
climate change differ markedly among species according 
to their adaptations to seed dispersal and growth forms, 
and that it is possible to predict the genetic consequences 
of range reduction by combining species modelling 
approaches with prior knowledge on species traits and/ 
or F S t estimates. 

2. MATERIAL AND METHODS 

(a) Species 

We analysed 27 plant species typically occurring in the bio- 
climatic zones at the tree line and beyond, i.e. the alpine and 
arctic zones. Information on the following species traits 
was compiled from the literature and databases: dispersal 
adaptation, growth form, pollination mode, breeding system, 
northernmost bioclimatic zone where it occurs, temperature 
tolerance and current geographical distribution type (see the 
electronic supplementary material) . 

(b) Loss of genetic diversity 

We analysed 73-958 individual plants from 14 to 131 local 
populations of each species and assessed levels of genetic 
diversity from AFLP data (78-334 markers per species, 
see the electronic supplementary material). For 24 of the 
species, full details of data collection and genetic structuring 
have been published elsewhere (electronic supplementary 
material, table SI). 

To estimate loss of genetic diversity expected as a conse- 
quence of range loss, we first divided the total study area 
according to an arbitrary grid of 500 x 500 km cells using 
ArcMap v. 9.2 and the Lambert azimuthal-equal area projec- 
tion of the Northern Hemisphere (figure 1). This was 
performed to account for differences in sampling intensity 
in different areas. For each species, only grid cells containing 
sampling localities were retained. The genetic consequences 
of range reduction were estimated as the loss of AFLP mar- 
kers occurring from randomly removing an increasing 
number of grid cells. This procedure was repeated 1000 
times to create a look-up table for minimum, maximum, 
mean and median number of markers lost for increasing 
numbers of grid cells removed. 

We expressed the number of AFLP markers lost relative to 
genetic diversity among individuals. Considering that a species 
would have lost all its genetic diversity when the last remaining 
population consists of identical individuals, we defined total 
loss of diversity as the total number of markers minus average 
number of markers per individual. The estimated loss of gen- 
etic diversity (Gi oss ) was thus expressed relative to this total loss 
as: G loss = average number of markers lost/(total number of 
markers — average number of markers per individual). We 
investigated the sensitivity of our estimates to variation in 
grid sizes and to geographical patterns in loss of genetic diver- 
sity (loss of adjacent grid cells after starting at a random point, 
see the electronic supplementary material) . 




Figure 1. Estimating loss of genetic diversity and range 
reduction, exemplified by data for Vaccinium uliginosum. 
(a) DNA sampling points, the 500 x 500 km grid overlaying 
the sampling points to adjust for variation in sampling inten- 
sity, and distribution sampling points used for modelling 
present and future distribution. GBIF, Global Biodiversity 
Information Facility. Dark grey shows current distribution. 
(6) Potential present and future (year 2080) distribution habi- 
tats overlaid to show lost (red), stable (purple) and future new 
habitat (blue). In this example, a 26% range reduction was 
estimated for the A2 emission scenario and CCM3 global cir- 
culation model, and there were 53 grid cells of 500 x 500 km 
that contained samples. The predicted loss of 26% of the 
range corresponded to a loss of 14 grid cells. 

(c) Range reduction and range gain 

To estimate range reduction and gain, we used SDMs relat- 
ing observed species occurrences to environmental variables 
[21]. These models are reported to be of moderate to good 
quality with respect to reconstructing recent and Holocene 
past distributions [22,23]. We selected three climatic variab- 
les regarded as important in determining plant distributions 
[24] : annual sum of precipitation, mean maximum tempera- 
ture of the warmest month and mean minimum temperature 
of the coldest month. These variables for current (1961- 
1990) and future (2071-2100) conditions were extracted 
from the Climate Research Unit (http://www.cru.uea.ac.uk/) 
and the Worldclim (http://www.worldclim.org/) data centre 
at a 10' resolution. Future climates were represented by 
Community Climate Model version 3 (CCM3) and Hadley 
Centre Coupled Model version 3 (HadCM3) GCMs for 
A2 (+3.2°C) and B2 (+1.0°C) emission scenarios [25]. 
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Observed distributions were digitized at the circumboreal/cir- 
cumarctic scale from Hulten & Fries' [26] distribution maps. 
Because these polygon maps represented the extent of occur- 
rence of the species, we selected a number of random 
presences inside these polygons (resulting in one point for 
1000 km 2 ). We also added validated occurrence points from 
the Global Biodiversity Information Facility (http://www.gbif. 
org/)., and from our own field sampling locations (electronic sup- 
plementary material, table SI). Selected climatic variables were 
extracted for these points. We created spatially random pseudo- 
absences (same number as presences) outside of these climatic 
boundaries [27] . All presence and absence locations were sum- 
marized in a 7.5 km resolution grid (North Pole Lambert 
azimuthal-equal area projection, 1.4 million cells) to decrease 
the pervasive effects of spatial autocorrelation. Current and 
future climate were then resampled in the same 7.5 km grid to 
run the projections. 

SDMs were calibrated within the biomod package in R [28], 
which capitalizes on seven widely used techniques (generalized 
linear models, generalized additive models, classification tree 
analysis, boosted regression trees, random forest, multiple 
adaptive regression splines and mixture discriminant analysis) 
to provide an ensemble of spatial projections. Model calibration 
was performed on a random sample of the data (70%) and 
model evaluation was carried out on the remaining 30 per 
cent with the true skill statistic (TSS) [29] and the area 
under the receiver-operating characteristic plot [30] . Calibrated 
models were then used to project current and future suitable 
climatic habitats over the Northern Hemisphere north of 
20°N latitude. We transformed the current occurrence prob- 
abilities into presences/absences using two methods; the 
thresholds that maximize (i) both the percentage of presences 
and absences correctly predicted, and (ii) the TSS. We used 
the same thresholds to convert future projections. To investigate 
the uncertainty coming from the overall modelling approach, 
we used the projections from the seven different SDMs, trans- 
formed into presence -absence using the two above-mentioned 
approaches (7x2=14 projections for current climate and 
each pair global circulation x emission scenario) to estimate 
percentage of range reduction, range gain and range change. 

Because the above-mentioned models can be overly com- 
plex, we also performed a simple rectilinear envelop model 
(surface range envelop [31]). The envelope is defined by 
identifying maximum and minimum values for each variable 
from observed presence of a species. Any presence with all 
variables falling between these maximum and minimum 
limits was included within the range that depicts the climatic 
conditions within which the species have been recorded. 
Being less constrained, the potential range is usually larger 
than those estimated by more complex models. 

(d) Loss of genetic diversity in the year 2080 

The proportion of genetic diversity likely to be lost in the year 
2080 was determined by evaluating the number of grid cells 
corresponding to the estimated range reduction (see figure 1), 
and then using the look-up table described earlier to deter- 
mine the expected loss of genetic diversity for this value. 
By choosing to use only predicted range loss, we ignore 
that the levels of genetic diversity in the new range may 
increase owing to increased mutation rate or introgression 
[32]. We considered this to be suitable given that genetic 
diversity in a new range has been shown to typically represent 
a subset of the genetic diversity in the source populations 
owing to bottlenecks and founder events [11,12,32], and 



given that we are considering a short time span. Based on a 
similar argument, we also ignore that genetic diversity may 
be retained if genotypes in extirpating populations manage 
to escape into the stable parts of the range; studies have 
shown that populations from southern portions of a geo- 
graphical range may contribute little to genetic diversity of 
more northern populations [33-35]. 

(e) Statistical analyses 

For each species, genetic differentiation among all sampled 
populations (F ST ) was estimated by an analysis of molecular 
variance using the software Arlequin v. 3.0 [36]. 

We tested for correlations among species' traits using 
linear models for correlations involving the continuous vari- 
ables like northernmost bioclimatic zone and temperature, 
and x 2 tests for 2 x 2 contingency tables to test the strength 
of association between categorical variables. We found a sig- 
nificant correlation between growth form and dispersal 
adaptation (electronic supplementary material, table S2). 
As dispersal has higher impact on genetic diversity than 
growth form [19,20], we kept dispersal adaptation. There 
was also a significant correlation between northernmost 
bioclimatic zones and both dispersal adaptations and breed- 
ing system (electronic supplementary material, table S2). 
As the thermal tolerance of the species was taken into 
account through maximum and minimum temperatures, we 
omitted the northernmost bioclimatic zone from further ana- 
lyses. There was a marginal significant interaction between 
breeding system and distribution, but this was ignored. 

Expected loss of genetic diversity, which was expressed as 
proportion, was arcsine (square root) transformed in all ana- 
lyses, which according to Shapiro -Wilk tests normalized 
distribution in all cases (p > 0.05) except for at 10 per cent 
range reduction. 

Linear models were used: (i) to investigate the effect of 
dispersal adaptations, pollination mode, breeding system, 
temperature tolerance and current geographical distribution 
type on loss of genetic diversity given a certain range loss, 
(ii) the usefulness of genetic differentiation (F ST ) as a predic- 
tor for the rate of genetic loss, and (iii) the effect of range 
reduction by 2080 on loss of genetic diversity. Model selec- 
tion was performed by Akaike's Information Criterion 
corrected (AICc) for small sample size ([37]; electronic sup- 
plementary material, table S3). This model selection 
criterion balances model complexity and model fit. The 
AAICc of the best-fit model is zero, and models with 
AAICc < 2 have some support [37]. Candidate models 
included models with single predictors and with all different 
combinations of two predictors (electronic supplementary 
material, table S3). The residuals of all selected models 
were checked using diagnostic plots to see if they satisfac- 
torily met the assumptions of linear models (stable 
variance). Some outliers were discovered, but the same or 
stronger correlations were found when we excluded the out- 
liers. To check for possible bias owing to unequal numbers of 
populations, individuals per populations or number of AFLP 
markers, each of these variables was added to the best candi- 
date model as additional explanation variables. All analyses 
mentioned above were performed in R v. 2.12 [38]. 

A possible problem with the above-mentioned analyses is 
that species are treated as statistically independent units. 
This approach compares the extant species situated at the 
tip of the phylogenies (TIP). As plant traits are likely to 
be correlated for phylogenetically related species, analyses 
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Figure 2. Estimated loss of genetic diversity as a function of decreasing range for 27 northern plant species. The bold line refers to 
the median; the dark grey shaded area refers to 50% CI; the light grey shaded area refers to 90% CI; and the dashed lines refer to 
minimum and maximum loss of genetic diversity. Vertical red lines show minimum and maximum range reduction expected by 
the year 2080 by any of seven species distribution models., two emission scenarios and two global circulation models (see §2). 



taking into account phylogeny may reveal different results. 
Therefore, we also tested the effect of traits using phyloge- 
netic-independent contrasts (PICs) [39]. We constructed a 
phylogenetic tree (electronic supplementary material, 
figure S2) using the online software Phylomatic (http:// 
www.phylodiversity.net/phylomatic/phylomatic.html) and 
choosing the angiosperm consensus tree [40]. A polytomy 
in Ericaceae was resolved according to Kron et al. [41]. 
The gymnosperm Juniperus communis was regarded as 
sister to all angiosperms. All branch lengths were assigned 
a value of 1. The categorical variables such as dispersal, 
range, breeding, pollination and life form were coded as 
dummy variables (0, 1). The PIC analyses were run in 
Compare v. 4.6b [42]. A rough 95% CI for the regression 
slopes was estimated as ±1.96s.e. Effects were considered 
significant, if the CI excluded zero. 



3. RESULTS 

In all cases, the relationship between loss of genetic 
diversity and loss of range was nonlinear, indicating 
that the majority of genetic markers are shared among 
geographical regions (figure 2). The rather narrow 90% 
CI observed for the majority of species indicate that 
loss of genetic diversity was in most cases rather indepen- 
dent on the order in which parts of the range were lost, 
and that loss could be estimated with high precision. 
For some species, however, uncertainties were high, 
especially Arabis alpina and Ranculus glacialis (figure 2). 
The median, minimum and maximum loss of gene- 
tic diversity remained largely similar in the analysis 
modelling loss of geographically adjacent areas, but the 
CI increased (see the electronic supplementary material, 
figure SI). 
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Loss of genetic diversity was best explained by disper- 
sal adaptation alone, or by dispersal adaptation combined 
with minimum temperature of the coldest month, or with 
total distribution range according to the AICc (electronic 
supplementary material, table S3). These results were 
robust to differences in number of populations, mean 
number of individuals per population and total number 
of polymorphic markers (see the electronic supplemen- 
tary material). Our estimates indicated that species 
without adaptations to long-distance dispersal (and/or 
herbaceous species) will lose genetic diversity at about 
twice the rate of species adapted to long-distance disper- 
sal by animals or wind (and/or woody species; figure 2 
and table 1). The results were similar for TIP and PIC 
analyses (table 1). However, dispersal adaptation was 
strongly correlated with growth form in our set of species 
(table 2 and electronic supplementary material, S2); thus, 
the effect of these two traits could not be distinguished. 
Importantly, but not unexpectedly, the rate of genetic 
loss with range reduction was estimated to be higher for 
species with higher genetic differentiation among popu- 
lations (F S t; table 1). Most of the species lacking 
adaptation to dispersal had F S t values greater than 0.5, 
whereas all but one species adapted to long-distance 
dispersal had F S t values less than 0.5 (table 2). 

The prediction power of the SDMs was generally high 
(area under the curve > 0.98, TSS > 0.85, electronic 
supplementary material, table S4) . For 1 8 of the species, 
the predicted range reduction (median based on seven 
techniques) exceeded 40 per cent for at least one emission 
scenario and circulation model (table 2) . Range reduction 
was on average higher under emission scenario A2 
('business as usual', 36-43% reduction) than under B2 
('reduced C0 2 emission', 26-32% reduction), as 
expected as the A2 scenario anticipates a more severe cli- 
mate change than the B2 scenario (table 2). The range 
gain was generally considerably lower than the range 
reduction, and on average the range change was —24 
per cent for A2 CCM3, -30 per cent for A2 HadCM3, 
-16 per cent for B2 CCM3 and -22 per cent for B2 
HadCM3 (electronic supplementary material, table S5). 

According to the median of seven techniques, all 
species except two very widespread, abundant and bird- 
dispersed ones (J. communis and Vaccinium uliginosum) 
were predicted to lose some of their present genetic diver- 
sity by 2080 (table 2). The estimated loss of genetic 
diversity in Micranthes foliolosa was 30 per cent, and six 
species were estimated to lose more than 20 per cent of 
their genetic diversity according to at least one scenario 
(table 2). However, there was a large gap between mini- 
mum and maximum estimated range reduction, and all 
species were expected to lose range under some models 
(electronic supplementary material, table S5). In the 
worst case scenario, assuming that the model estimating 
the maximum range reduction will be realized, and that 
the corresponding loss of genetic diversity will be at the 
maximum value (crossing point between right red bar 
and upper dashed line in figure 2), all species were 
expected to lose some genetic diversity, one-third of 
them greater than 50 per cent. 

The median estimated loss of genetic diversity differed 
only slightly among the emission scenarios and GCMs 
(table 2). Species expected to lose more of their range 
were also prone to more severe genetic loss: range reduction 



alone explained 66-74% of the variation in estimated loss 
of genetic diversity among species (not shown) . 

4. DISCUSSION 

Our results show that species traits which are well known 
to influence patterns of genetic diversity within and 
among populations, also affect the predicted loss of gen- 
etic diversity, and thereby the susceptibility of species to 
genetic depauperation under range reduction. As our 
modelling predicts range reduction and loss of genetic 
diversity by 2080 in all 27 species according to at least 
one scenario, it can be expected that climate warming 
will have a major impact on the future range sizes and 
levels of genetic diversity in northern plant species. 
Because genetic diversity is important not only for 
species' persistence and evolutionary potential [3,4], but 
also for community structure and ecosystem resilience 
[43], climate change-induced loss of genetic diversity 
may be expected to affect all levels of biodiversity. 

(a) Vulnerability to loss of genetic diversity 
The precision of our predicted loss of genetic diversity lar- 
gely depended on the geographical distribution of genetic 
diversity within the species. For species with a rather even 
level of genetic diversity throughout most of its distribution 
range, as e.g. in Carex bigelowii and S. herbacea [18,44,45], 
the prediction intervals are fairly small with the minimum 
and maximum predicted loss of genetic diversity close to 
the median, indicating that the effect of range reduction 
is rather independent on which part of the range is lost 
(figure 2). However, for species such as those with high 
levels of genetic diversity in southern alpine areas and 
hardly any diversity in northern areas, as observed in 
A. alpina and R. glacialis [46,47], the future loss of genetic 
diversity will strongly depend on which part of the range is 
lost. As southern populations are most likely to get lost 
(cf. figure 1), the expected loss of genetic diversity is 
likely to be at the maximum estimates (figure 2), and 
thus cause severe genetic depauperation. 

Dispersal adaptations appear to be important in deter- 
mining the rate of loss of genetic diversity. This is not 
surprising, as dispersal adaptations influence the level of gen- 
etic differentiation among plant populations [19,20,48]. 
Also among our species, those adapted to long-distance dis- 
persal had lower F S t values than the species lacking such 
adaptations. However, given the strong correlation between 
dispersal adaptations and growth form in our dataset, the 
differences we observed among short- and long-distance 
dispersers may be based on a combination of these two 
species traits. Growth form may indeed also contribute to 
shaping genetic patterns, as woody species (especially 
trees) often have a larger stature and lower population 
densities, which could result in higher pollen and seed dis- 
persal [19]. A recent meta-analysis showed, however, that 
woody and herbaceous species may experience similar loss 
of genetic diversity, a consequence of habitat fragmentation 
[14]. Woody species may also suffer more from inbreeding 
depression than herbaceous plants, and thus have a higher 
selection pressure than inbreds [49]. In addition, many 
species among our long-distance dispersers occur abundan- 
tly as vegetation dominants with very high seed outputs, 
which may further facilitate long-distant gene flow. The 
typically high levels of gene flow among populations of 
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Table 2. Estimated range reduction CR re d.3%3 median of seven techniques) and associated estimated loss of genetic diversity 
(Glossy %) for 27 northern plant species., according to the A2 and B2 emission scenarios and the CCM3 and HadCM3 global 
circulation models. (Dispersal: S, short-distance; L, long-distance. Growth form: H, herbaceous; W, woody. Genetic 
differentiation among populations (Fst)-) 
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long-distance-dispersed woody species [19,20] create a gen- 
etic pattern which make species less vulnerable to loss of 
genetic diversity during range reduction. As loss of genetic 
diversity is expected to adversely affect the ability of popu- 
lations to evolve and cope with environmental change and 
thereby increase the risk of extinction [50,51], short- 
distance-dispersed herbs may be expected to enter the 
extinction vortex more rapidly. Using knowledge on species 
traits may thus help us forecast which species are at risk. 

The strong correlation we found between F S t and loss 
of genetic diversity, suggests that F S t values can be used 
as a proxy for predicting the genetic vulnerability of 
species to climate change. However, as estimates of F ST) 
G S t or related measures of population differentiation 
depend on the type of genetic markers used, the geogra- 
phical scale of the study, the number of individuals and 
populations sampled, and the estimator used [52-54], 
the actual F S t estimates should be evaluated in the light 
of these sources of bias. 

Our AFLP-based predicted loss mainly represents 
loss of neutral genetic diversity [55], whereas species' 
abilities to adapt and survive climate change mainly 
depend on adaptive genetic variation. Although methods 



for estimating range-wide levels of adaptive genetic vari- 
ation are under development [56], it is still not feasible 
to apply them to approximately 10 000 samples of 27 
different species as studied here. Neutral genetic variation 
can be expected to be lost at a rate similar to that of adap- 
tive genetic variation [57], but the differentiation among 
populations is typically higher when based on adaptive 
than on neutral genetic diversity [58,59]. It is therefore 
likely that our estimates are underestimating the actual 
loss of evolutionary potential. 



(b) Prediction for 2080 

Overall, the range changes modelled for these northern 
species in 2080 (electronic supplementary material, table 
S5) were slightly higher than those modelled for temperate 
tree species under the HadCM3 GCM and A2 (30% versus 
22%) or B2 (22% versus 19%) emission scenarios [60], 
and that modelled for 84 Danish temperate species (29% 
and 12% for scenarios A2 and B2, respectively) [61]. In 
addition, alpine species are expected to lose more habitat 
than boreal species [10]. Thus, it is reasonable to suggest 
that susceptibility to climate warming-induced range 
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reduction and associated loss of genetic diversity will be 
highest in species restricted to cold climates. 

Currently, The International Union for Conservation of 
Nature [62] red list criteria do not take genetic diversity 
into account, but species with a population reduction 
expected in the future (up to a maximum of 100 years) 
should be listed as critically endangered, endangered or 
vulnerable if the expected range reduction is greater than 
or equal to 80 per cent, greater than or equal to 50 per 
cent, or greater than or equal to 30 per cent, respectively. 
As many as 26 of the 27 rather common and widely distrib- 
uted species studied here fall within these categories 
according to at least one emission scenario, modelling tech- 
nique and global distribution model (figure 2). Given the 
high number of species that are expected to become vulner- 
able or threatened by 2080 [10], it will become difficult to 
prioritize species for conservation. As our study indicates 
that the genetic consequences of range reductions are 
remarkably different among species, genetic parameters 
should be considered in future management assessment. 

SDMs are widely used for management purposes [63]. 
The prediction power of these models in our dataset was 
high, and we show that they may also, combined with 
genetic data, be a useful tool for estimating expected loss 
of genetic diversity. Nevertheless, these models do not 
take into account mechanisms such as changes in biotic 
interactions, potential for rapid adaptations or time lag. 
Although some newly developed spatially explicit models 
are able to account for these factors [64,65], they are not 
ready to be implemented for a large set of species and to 
be run over the whole arctic- alpine region. When species 
shift ranges, some populations may persist at the rear edge 
and conserve local genetic variants. If this is the case, then 
our approach may overestimate habitat loss and associated 
loss of genetic diversity. Research on rear-edge populations 
is still limited [34], and their response to climate change 
may be challenging to model [2]. Populations may persist 
if they are able to colonize adjacent microhabitats in hetero- 
geneous landscapes [2,34], if they are able to endure long 
periods without recruitment (owing to, e.g. clonal growth, 
persistent seed bank, long lifespan, etc.) [2,34], and/or 
because they adapt to the changing climate [7,66,67]. We 
assume that switching of microhabitats is equally likely in 
short- and long-distance-dispersed species. Extinction lags 
are more pronounced in long-lived woody species than in 
herbaceous species [9], but this might be outweighed by 
the faster evolutionary rate in herbaceous than woody 
species [68]. Also, as good dispersal ability is likely to 
enhance both the ability of range change and adaptive evol- 
ution [67], we think that any uncertainties in the models are 
likely to increase rather than decrease the differences we 
observed in susceptibility to loss of genetic diversity in 
woody, long-distance-dispersed compared with herbaceous, 
short-distance-dispersed species . 

5. CONCLUSIONS 

Our study shows that F S t values, which are available for a 
wide range of species, as well as species traits such as dis- 
persal adaptation and growth form, can be used to predict 
a species' susceptibility to loss of genetic diversity follow- 
ing climate change. As it is important to assess and 
commence management actions before genetic diversity 
is lost [69], we advocate to combine SDM with data on 



genetic differentiation and/or species traits to predict 
which species are at highest risk of losing genetic diversity 
in a changing climate. Such an approach will facilitate 
rational prioritization of conservation efforts. 
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